Scaling of Berry's Phase Close to the Dicke Quantum Phase Transition 
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We discuss the thermodynamic and finite size scaling properties of the geometric phase in the 
adiabatic Dicke model, describing the super-radiant phase transition for an A'^ qubit register coupled 
to a slow oscillator mode. We show that, in the thermodynamic limit, a non zero Berry phase is 
obtained only if a path in parameter space is followed that encircles the critical point. Furthermore, 
we investigate the precursors of this critical behavior for a system with finite size and obtain the 
leading order in the 1/N expansion of the Berry phase and its critical exponent. 
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A considerable understanding of the formal description 
of quantum mechanics has been achieved after Berry's 
discovery l"! of a geometric feature related to the motion 
of a quantum system. He showed that the wave function 
of a quantum system retains a memory of its evolution 
in its complex phase, which only depends on the "ge- 
ometry" of the path traversed by the system. Known as 
the geometric phase, this contribution originates from the 
very heart of the structure of quantum mechanics. His- 
torically, the first implicit derivation of geometric phase 
is to be found in the work of H. C. Longuet-Higgins et 
al. [31 and G. Herzberg and H. C. Longuet-Higgins 3] on 
the study of the E ^ e Jahn- Teller problem in molecular 
systems. They showed that the Born-Oppenheimer wave 
function, if required to be real and smoothly varying, un- 
dergoes a sign change (a tt geometric phase shift) when 
the nuclear configuration is brought around a crossing 
point of the energy levels, i.e. a point where the en- 
ergy levels become degenerate. The widespread interest 
in geometric phases is unquestionably due to the work of 
Berry 'l'| who put the concept of geometric phase in an 
unexpectedly general scenario, showing that it is a prop- 
erty of the evolution of any quantum mechanical system. 
It is then not too surprising how this idea has been so 
broadly studied and applied to a variety of contexts |j,|5|. 

An apparently unrelated area in which a connection 
with Berry Phase (BP) has been recently drawn is the 
study of Quantum Phase Transitions (QPT) in many- 
body systems [a, LD, IS|- It is a well known fact that 
QPT |9| are accompanied by a qualitative change in the 
nature of correlations in the ground state of a quantum 
system, and describing these changes is clearly one of the 
major interests in condensed matter physics. In partic- 
ular, critical points are associated with a non-analytical 
behavior of the system energy density [Sj . It is expected 
that such drastic changes are reflected in the geometry of 
the Hilbert space. The geometric phase is able to capture 
singular behaviors of the wave function, and is therefore 
expected to signal the presence of QPT. The appearance 
of non-trivial BP in presence of criticality is also heuris- 



tically motivated by the fact that QPT arise in corre- 
spondence of energy level crossings (in the thermody- 
namical limit) between ground and excited states 9]. As 
implicitly suggested by the early work of H. C. Longuet- 
Higgins et al. I2] , level crossings generate singularities in 
the curvature of the Hilbert space. Having the ability 
to adiabatically change the parameters of the system so 
as to encircle these singularities, can reveal the critical 
dependence of geometric phases 6,7]. Indeed, due to its 
topological character, the BP provides a tool to detect 
and probe criticalities without taking the system through 
a phase transition. This has been recently considered in 
ref. |6|, where a (1+1) XY spin model is shown to display 
a topological behavior of the BP when a path that circu- 
lates (or not) a region of XX criticality is taken. Zhu have 
analyzed further the connection between BP and QPT, 
obtaining a scaling relation for the geometric phase J3]- 
Very recently, Hamma has demonstrated the topological 
character of the BP across regions of QPT for a general 
model of many body system 7] . 

These results have been achieved in the thermody- 
namic limit and it would be interesting, in this context, 
to understand whether and how precursors of QPT for 
finite size systems appear in the geometric phase. In this 
letter, we analyze this problem for the case of the critical 
behavior of the Dicke Model (DM) and investigate the 
topological and scaling properties of the geometric phase 
across a region of criticality. We consider a system which 
consists of N two-level systems (a qubit register or an 
ensamble of indistinguishable atoms) coupled to a single 
oscillator (bosonic) mode. The Hamiltonian is given by 
(in unit such that h = c = 1) 



H = uja'^a + ASa; + —^{a'^ + a)Sz 

\/N 



(1) 



where A is the transition frequency of the qubit, oj is the 
frequency of the oscillator and A is the coupling strength. 
The qubit operators are expressed in terms of total spin 
components Sk — J2i=i '^L where the cr^ 's (fc — x, y, z) 
are the Pauli matrices used to describe the i-th qubit. A 



7r/2 rotation around the y axis shows that H is canoni- 
cally equivalent to the Dicke Hamiltonian [lO|, including 
counter-rotating terms. 

After the first derivation due to Hepp and Lieb [llj , the 
thermodynamic properties of the DM have been studied 
by many authors Hi E H E E [13 . In the thermo- 
dynamic limit {N ^ oo), the system exhibits a second- 
order phase transition at the critical point Ac — yKuj/2, 
where the ground state changes from a normal to a super- 
radiant phase in which both the field occupation and the 
spin magnetization acquire macroscopic values. The con- 
tinued interest in DM stems from the fact that it dis- 
plays a rich dynamics where manv non-classical effects 
have been predicted [l^, lid . 12(1 I21I l22| . and from its 
broad range of applications [22|. Investigations of the 
ground state entan glem ent of the DM have been recently 
performed |2J, |25|, |20| , pointing out a scaling behavior 
around the critical point. 

The aim of present work is to investigate the DM in 
the adiabatic regime (A ^ uj) and to demonstrate the 
topological character of the BP for this case. Further- 
more, we show that the geometric phase obeys scaling 
relation at the critical point for a system with finite size 

In order to generate a Berry phase we change the 
Hamiltonian by means of the unitary transformation: 



U{(j)) = exp 



2 



(2) 



As the qubits are indistinguishable, it is easy to prove 
that the ground state can be expressed as a direct prod- 
uct of N identical factors, 

lx((Z, 0)) = lx(g, ^))i ® lx(9, 4>))2 • • • ® |x(g, c^))n (6) 
Each component can be written as 

|x((?,0)).=sin^|r).-cos^e<U)„ (7) 

where | t)i and | |)i are the eigenstates of af with eigen- 
values ±1, and where 
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ND 



(8) 
(9) 



Here, E{q) is related to the energy eigenvalue of Eq. Q 
as 



Ei{q)^-NE{q) = -N\D 



L^q^ 



(10) 



In the Born-Oppcnheimer approach, this energy eigen- 
value constitutes an effective adiabatic potential felt by 
the oscillator together with the original square term: 



Viiq)^^[q'-NEiq)] 



(11) 



where is a slowly varying parameter, changing from 
to27i-. 

The transformed Hamiltonian can be written as 



H(0) = C/t (0)i7[/(0) = -[p^ + q^ + G-S] (3) 



where G 



is an effective vec- 



tor field. Here, D — 2A/ijj and L = 2^/2X/lj are di- 
mensionless parameters and the Hamiltonian of the free 
oscillator field is expressed in terms of canonical vari- 
ables q = {a^ + a)/\/2 and p = i{a^ — a)/\/2 that obey 
the quantization condition [q,p] = i- 

In the adiabatic limit [23, I23 , where we assume a slow 
oscillator and work in the regime D ^ 1, the Born- 
Oppenheimer approximation can be employed to write 
the ground state of H{(p) as: 



li^t. 



dq ip{q)\q) (g> |x(q, 0)) 



(4) 



Here, |x(q, 0)) is the state of the "fast component"; 
namely, the lowest eigenstate of the "adiabatic" equation 
for the qubit part, displaying a parametric dependence on 
the slow oscillator variable q, 



G-S\x{q,^))^Ei{q)\xiq,^)), 



(5) 



Introducing the dimensionless parameter a = L^/2D, 
one can show that for a < 1, the potential Vi{q) can be 
viewed as a broadened harmonic well with minimum at 
q = and V;(0) = — A^A. For a > 1, on the other hand, 
the coupling with the qubit splits the oscillator poten- 
tial producing a symmetric double well with minima at 
±q^ = ± V^V^^^T with Viiq^) = -i^ (a -f i). 

As the last step in the Born-Oppenheimer procedure, 
we need to evaluate the ground state wave function for 
the oscillator, (j)o{q), that has to be inserted in Eq. Q 
to obtain the ground state of the composite system. 
This wave function is the normalized solution of the one- 
dimensional time independent Schrodinger equation 



2d^ 



Had^poiq) = - 



+ Viiq)]M<l)^£oM<l), (12) 



where £0 is the lowest eigenvalues of the adiabatic Hamil- 
tonian defined by the first equality. 

Once this procedure is carried out for every value of the 
rotation angle (p, the BP of the ground state is obtained 
as 



7 = i d(f>{ip, 



^t. 



l'2-rr 

dqVoiq) / d(l)A{q,, 
Jo 



(13) 
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FIG. 1: Numerical results for the scaled BP as a function of 
the parameter a, for D = 10 and for different values of A^, 
in comparison with the result for N ^ oo. The BP increases 
with the coupling, and, in the thermodynamic limit, it dis- 
plays a cusp at the critical value a = 1. The inset shows the 
derivative of the BP with respect to a. 



where the ^-parametrized connection is 
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d(l 

Lq_ 

Vn 



ND 



dC „„„2 /3 



xii'i,(t>) ) = -^-n^'^^ 



COS( 



2E{q) E{q) - ^ COS 



2 

(14) 



Substituting this expression into Eq. I|13|l , straightforward 
calculations lead to 



7 = A^TT 14- 



N 



where 



N 



'pUq) 



D 



dq , 



In the thermodynamic limit, one can show 



N 



-1 (a<l) 



and thus, for N ~* 00, the BP is given by 



1_ 
N 



N- 



(a < 1) 

7r(l-i) (a>l) 



(15) 

(16) 
that 

(17) 

(18) 



Numerical results for the scaled BP are plotted in Fig. (^ 
as a function of the parameter a, for Z? = 10 and for dif- 
ferent values of N, in comparison with the result for the 
thermodynamic limit. One can see that the BP in- 
creases with increasing the coupling strength and, in the 
thermodynamic limit, its derivative becomes discontinu- 
ous at the critical value a = Oc = 1. This results are in 
accordance with the expected behavior of the BP across 
the critical point. Notice that, in the thermodynamic 
limit, a non-trivial Berry phase is only obtained when 
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FIG. 2; A qualitative illustration of the paths followed by the 
parameters of the Hamiltonian due to the application of U{<j)). 
The paraboloid corresponds to the value a = L /2D = 1, 
for which the Hamiltonian shows a critical behavior. If the 
parameters follow a path, e.g. Ci, encircling the paraboloid, 
then the system acquires a non-trivial BP, which tends to ir 
for a ^ 1. As seen from figure Q path C2 gives rise to a zero 
Berry phase (in the thermodynamical limit). 



a region of criticality is encircled, as for the path C'l in 
Fig. |2] Indeed, in the enlarged parameter space gener- 
ated by the application of the unitary operator U{(j)) of 
Eq. |(2Jl, the critical point corresponds to the paraboloid 
a — ^ ~ I. As the radius of the path is determined 
by a, one can see that, in the limit N -^ 00, the BP is 
zero in the normal phase (a < 1) and is non-zero in the 
super-radiant phase, i.e. if the path encloses the critical 
region. 

We now investigate the scaling behavior of BP at the 
critical point by a finite size scaling approach. In order to 
obtain an analytic estimation of BP as a function of N, 
we expand the adiabatic potential in Ea. (|12|) and obtain 
an anharmonic oscillator potential 



Ui{q) = -Vi{q) ~ -ND + (1 - a)q' + 

UJ 



2ND 



(19) 



The eigenvalue problem defined by this potential can be 
solved with the help of Symanzik scaling |28l |22j . This 
is done by rewriting Ea. H12|l into the equivalent form 






ipo{x;C)^eoiC)M^-X) (20) 



where the scaled position is a; = g I „° „ j 
,2/3 



1/6 



while ( 



{^^^Y'"^ (ac — a). Finally, the energies in Ea. l(T^ and 
(I20II obey the scaling relation 



-SQ{a,ND) = -ND 



f a^ 



\2ND 



1/3 



eo(C). (21) 



Since C ^ at the critical point, we can consider the 
x^ term to be a perturbation and employ the Rayleigh- 
Schrodinger perturbation theory. This yields the expan- 
sion eo(C) = Sn=o CnC"? where the coefficients c„ can be 



log 




-10- 



log2(N) 



FIG. 3: Scaling of the BP as a function of A^ at the crit- 
ical point a = 1, for D = 10. For easy of comparison, the 
continuous plot shows the analytic expression of Eq. 12311 . 



argument motivates the need to explore the use of BP as 
a tool to signal and investigate critical features of certain 
models. However, strictly speaking, singularities in the 
energy density of many body systems only appear in the 
thermodynamic limit. It is therefore not obvious that in 
a finite scale regime, such a connection between BP and 
QPT can still be drawn. Studying the BP in this regime 
has clearly theoretical interest and obvious experimental 
motivations. In the case of the Dicke model, we have 
found that the geometric phase shows, at finite sizes, a 
precursor of the topological character which appears in 
the thermodynamic limit. Moreover, studying the scal- 
ing of the BP as a function of the system size, we have 
identified its critical exponent. 



obtained after solving the equation for a purely quartic 
oscillator. It is easy to get cq = eo(0) ~ 1.06036 and 
ci = /_°1 g' '/'§(?; 0)dq = e[,(0) ~ 0.36203. 

It can be shown that the coefficients of this expansion 
completely determine the average value of every physi- 
cal observable at the critical point l28|. In particular, a 
similar expansion applied to (Sx), allows us to write 



N 



-1 



2ci 



2co 



(2iVL')2/3 {2ND)'^/^' 



(22) 



Thus, we obtain the leading orders in the finite size scal- 
ing of the Berry phase as 
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N 



2ci 



2co 



_{2ND)y^ {2ND)^/\ 



(23) 



This expression shows how the scaled geometric phase 
goes to zero as A'^ increases and how the singular ther- 
modynamic behavior is approached at a — ac = I. This 
analytical result, and in particular the leading critical be- 
havior 'y/N ^ A^~2/3 jg confirmed in Fig. by compar- 
ison with the BP obtained numerically from Eqs. H15|) - 
(|16|l . In fact, including also the second order correction 
scaling as N^^^^, we are able to reproduce the numerical 
result even for small values of N. 

Besides the scaling relation at the critical point a = 
1, we can also obtain the leading 1/N correction to the 
thermodynamic limit of -^ for small and large values of 
a. Using the fact that the oscillator localizes around 
q = for a <C 1, while its wave function is split in two 
components peaked around ±qm for a 3> 1, we get 
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N 



1_ 
N 



N- 



2ND 
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WDo? 



(a<l) 
(a>l) 



(24) 



To conclude, we have shown that the behavior of the 
geometric phase in correspondence of the critical region 
for the Dicke Model, confirms an expected connection 
between BP and QPT. Indeed, BP and QPT share the 
common feature of both appearing in presence of a singu- 
larity in the energy density of the system. This heuristic 
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